SPACE-TIME FRACTIONAL EQUATIONS AND THE RELATED 
STABLE PROCESSES AT RANDOM TIME 
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Abstract. In this paper we consider the general space-time fractional equa- 
tion Yl'jLi>^j-§p^w{^ir ■ ■ ,^n;t) = -c^ {-A)!^ w{xi,- ■ ■ ,Xn;t), for Uj e 
(0,1], 13 e (0,1] with initial condition w{xi,--- ,x„;0) = YYj^iSixj). We 
show that the solution of the Cauchy problem above coincides with the distri- 
bution of the n-dimensional vector process S^^ {c?L'^^'"' (t)) , t > 0, where 
S^^ is an isotropic stable process independent from X"!' " '""'{t) which is the 
inverse of jW"!. '"^ {t) = T,T=i ^j"' H^i (t),t> 0, with H^i (t) independent, 
positively-skewed stable r.v.'s of order Vj. The problem considered includes the 
fractional telegraph equation as a special case as well as the governing equation 
of stable processes. The composition S^^ (c^i'^i' " (t)) , t > 0, supplies a 
probabilistic representation for the solutions of the fractional equations above 
and coincides for /3 = 1 with the n-dimensional Brownian motion at the ran- 
dom time £''1 ' (i), t > 0. The iterated process S,^^' " (t), t > 0, inverse 

to fi^i'-'^-Ct) = T.T^i^Y'^' iH'^^iiH''^ izH'^^ {■■■ rH'^Ht) ■■■))), t > 0, 
permits us to construct the process [c^Z"^' '"'"(*)), t > 0, the distribu- 
tion of which solves a space-fractional equation of the form of the generalized 
fractional telegraph equation. For r — > oo and (3 = 1 we obtain a distribu- 
tion, independent from t, which represents the multidimensional generalisation 
of the Gauss-Laplace law and solves the equation AjUi(a;i, ■ ■ • ,x„) = 

n?i^'''(^i> ■ ' ■ y^n)- Our analysis represents a general framework of 

the interplay between fractional differential equations and composition of pro- 
cesses of which the iterated Brownian motion is a very particular case. 



1. Introduction and preliminaries 

1.1. Introduction. The study of tlie relationships between fractional differential 
equations and stochastic processes has gained considerable popularity during the 
past three decades. In pioneering works simple time-fractional diffusion equations 
have been considered (see for example Fujita [8]) and its connection with stable 
processes has been established (see Orsingher and Beghin [13]; the reader can also 
consult Zolotarev [18] for details on stable laws). In such papers the authors have 
shown that the compositions of processes have distributions satisfying fractional 
equations of different form. The iterated Brownian motion Bi {\B2(t)\), t > 0, (with 
Bi and B2 independent Brownian motions) has distribution solving the fractional 
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equation (see AUouba and Zheng [1]) 

'{x,t) = -^—u{x,t), x€R,t>0. (1.1) 
dt2 25 dx^ 

as well as the fourth-order equation (see DeBlassie [5]) 

It has been shown by different authors (see Benachour et al. [3]) that the solution 
to the biquadratic heat-equation 

d 1 

-u{x,t) = --—^u{x,t), a;eM,t>0, (1.3) 

coincides with 



u{x, = E < , cos -r-m ^ ) (1-4) 

^ ' \ V2^|i3(t)| \2\B{t)\-l^ ( ^ ^ 

and appears as the distribution of the composition of the Fresnel pseudoprocess 
with an independent Brownian motion (see Orsingher and D'Ovidio [14]). 
When the fractional telegraph equation 

' • 2X—]uix,t) ^ c^—u{x,t), xeR,t>0, (1.5) 



\dt^'' df" J ^ ' ' dx^ 

for 1/ e (0, 1], A > 0, c > 0, is considered, the solution of problem (1.5) ior u = ^ 
has been proved to coincide with the distribution of T{\B{t)\), t > 0, where T{t), 
i > 0, is a telegraph process independent from the Brownian motion B{t), t > 0. 
From the analytical point of view, equations similar to (1.5) have been studied in 
the form 

Qa qP ^7 

—u{x,t) + a-^u{x,t) = c^ — u{x,t)+(,'^u{x,t) + ip{x,t), xeR,t>0, (1.6) 

for a € [0, 1], /3 e [0,1], by Saxena et al. [15]. These authors have provided the 
Fourier transform of solutions of fractional equations of the form 

ai^^"(2^,i)H l-Q«+i g^e,„+i ^(^'^) ^ -g^u{x,t) + S,'^u{x,t) + Lp{x,t) (1.7) 

for ai, • • • , an+i G (0, 1) and > 0, (see Saxena et al. [16]) in terms of generalized 
Mittag-Leffler functions (but no probabilistc interpretation has been given to these 
solutions). Telegraph equations emerge in electrodynamics, in the study of damped 
vibrations, in the analysis of the telegraph process. Its multidimensional version 
appears in studying vibrations of membranes and other structures subject to fric- 
tion. Equations with many fractional derivatives emerge in the study of anomalous 
diffusions as pointed out by Saxena et al. [15, 16]. 

The symmetric stable laws have distribution satisfying the space- fractional equa- 
tion 

^u{x,t) ^ ~^u{x,t), xeR,t>0, (1.8) 
ot o\x\'^ 

where is the Riesz fractional derivative. For asymmetric stable laws the con- 
nection with fractional equations has been established by Feller [7] . The connection 
between fractional telegraph equations and stable laws has been established in a 
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recent paper by D'Ovidio et al. [G], in which the authors considered the multidi- 
mensional space-fractional extension of (1.5) 

^ + 2A^^ uix, t) = (-A)^ u{x, t), xe W\ t > 0, (1.9) 

for v £ (O, /3 e (0,1]. The solution to (1.9) subject to the initial condition 
u{x,0) — S{x) is given by the law of the composition of the form S^^ {c^L'^{t))^ 
t > 0, where Sl^{t), t > 0, is a n-dimensional isotropic stable vector process and 

L^it) ^ inf |s : Hf^is) + {2\)^ H'i{s) > t} (1.10) 

where Hf'^{t) and IH^it), i > 0, are independent positively-skewed stable processes, 
with V S (O, ^] . For (i — 1 the composition above takes the form of a Brownian 
motion at the delayed time L'^{t), t > 0. For = and n = 1 this establishes the 
fine distributional relationship 

T{\Bm'^^ B^c'dit)), t>0, (1.11) 

see D'Ovidio et al. [6]. 

In the present paper we consider the further generalization of the space-time 
fractional equation with an arbirtrary number of time-fractional derivatives 



Er=iA,w<,-,.™(a;,i) = -cM-A)'' (a;, t), xeW\t>0, 



if^i,-,i'„(a;,0) = 5{x), 



(1.12) 



c 



for Vj e (0, 1], /3 e (0, 1], Xj > 0, j = 1, • • • ,TO. The symbol stands for the 
Dzerbayshan-Caputo fractional derivative which is defined as 

^— /(t) = — ^ , — ds, m-Kiy <m,meN, (1.13) 

dt" ^ r (m - 1/) 7o {t - 5)^^+1-™ y ^ \ / 

for a function f E (M) (for fractional calculus the reader can consult Kilbas et al. 
[10]). The fractional Laplacian (— A)*^, /3 € (0, 1) is defined and explored in Section 
1.2.2 below. We show that the solution to (1.12) is given by the law of the process 

S^f (c^X^'i' - '""^(t)), t > 0, where 

L^^'-'^-^it) = inf {s > : ^''^•■■■'^'"(s) > t} (1.14) 

and 

m I 

^-i.-.-"(t) = ^Aj^ilf (i), t>0, (1.15) 
i=i 

for Hj^ ^ j — l^ - ■ ■ , TO, totally positively-skewed stable processes (stable subordina- 
tors) , of order i^j . In other words we show that the solution of a general space-time 
fractional equation (which includes reaction-diffusion equations, telegraph equa- 
tions, diffusion equations as very special cases) coincides with the distribution of 
a stable vector process taken at a random time L'^{t), t > 0, constructed as the 
inverse of the combination of independent stable subordinators. For the classi- 
cal Laplacian (/3 = 1) we have that the solution to (1.12) is the distribution of a 
Brownian motion at time L'^{t), t > Q. 
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We also prove that the law of the processes (1.14) and (1-15) are also solutions 
of fractional differential equations. In particular we show that 

= ^ (1.16) 

is the solution to the space-fractional problem for Vj e (0, 1) 

= EjLi ^i^^i'i, -,''n.(2^.*)> x>0,t>0 
/i,,,...,,„(rr,0) = <5(cc), (1.17) 
^,„(0,i) = 0, 

while the law of L"^'"' '"-^{t) solves 



(1.18) 



for Vj e (0, 1). In (1.17) and (1.18) the fractional derivatives must be meant in the 
Riemann-Liouville sense that is, for a function f E (M), 

(1;19) 

A section is devoted to the case of the fractional equation with two time deriva- 
tives of order a G (0, 1] and ly € (0, 1] with a ^ v, 

U^ + 2\y^) w'l^ {x, t) = (- A)^ {x,t), xeR-,t>0, 
\<,(a;,0) = 6{x), 

(1.20) 

which takes a telegraph-type structure for a = kv, k E N, kv < 1. The Fourier- 
Laplace transform of the solution of (1.20) for a ^ kv reads 

Jo JK" /.t*-'^ -I- 2A^'^ + cJ mil 

For fc = 2, 7i = l,/3 = l, we have the classical fractional telegraph equation studied 
in Orsingher and Beghin [12]. The Fourier transform of W2i/^iy (x, t) reads 

dxe^^'^W2,y^i,{x,t) = 



(1.22) 

where 771 and 772 are the solutions to + 2A^^ + c^^^ = and 

00 ^ 

^^,^(^)-E r(V.fc-fz?) ' ^,^>0,^GK, (1.23) 

is the two-parameter Mittag-Lefher function. For ly = 1, (1.22) coincides with the 
characteristic function of the telegraph process. For k — 3 and ^ < | in (1.21) we 
obtain explicitely the Fourier transform of the solutions in terms of Mittag-Leffler 
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functions and the Cardano roots A, B and C of the third order algebraic equations 
/^^'^ + 2X^1" + c2 ll^ll^'^ =0. For /c> 3 we can write 

(1-24) 

but the exphct evaluation of Z,; is, in general, impossible. 

In Orsingher and Beghin [13] n-times iterated Brownian motion 

3n{t) = Bi (1^2 (1^3 • • • (|B„+i(i)|) • • -1)1) , i > 0, (1.25) 

is considered and its connection with the fractional diffusion equation 

^uix,t) = 2^-^^uix,t) (1.26) 

investigated. Here we consider first the n-times iterated positively-skewed stable 
process jH'^^ with weights Xj > 0, j — 1, ■ ■ ■ , m, 

(i) = ^ iH''^- ( ^H-'^ ( ai?"^^ (• • • rH''^ {t) ■ ■ •))) , ^ > 0, (1.27) 

We construct the inverse of the process (1.27) as follows 

£;:^'-''^'"(i) = inf{s>0:io5:i'-'''"(s)>t}, t>0. (1.28) 

We show that the distribution of the composition 

B„(c2£,';i'--'^'"(t)), t>0, (1.29) 

where B„ represents the n-dimensional Brownian motion, is the solution to the 
Cauchy problem for i^j S (0, 1], re N, 

fE;iiA,^lre;:...,.,„(a^,t) = c2Ame;:...,.,„(a.,i), xeW\t>0, 
\n>^A,..^,^{x,0) ^ S{x). 
We show that for the number r of iterations tending to infinity 

law 

B„(c2£::i —-(*)) =^ Xra,n, (1-31) 

r — >■ oo 

where Xm.n is a r.v. independent from t and possesses density equal to 



dx (27r)t I c I II II ^1 ^ II I 

(1.32) 

where K^{x) is the modified Bessel function. For n = 1 the distribution (1.32) 
becomes the Gauss-Laplace law 



«.„(.)== ^%^e-^N. (1.33) 
2c 

Result (1.33) was obtained also in Orsingher and Beghin [13] and by a different 
approach for Ai = 1, Xj = for j > 2, c = |, was derived by Turban [17] as the 
limit of iterated random walks. 
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1.2. Preliminaries. 

1.2.1. One dimensional stable laws. Let us consider a stable process, say S'^(t), 
t > 0, ly E {0, 2], ly 1, for which, in general, 

Eg»€S""(t) _ g-(T|5|''t(l-i9sign(5)tan ig^) (1-34) 

where G [—1,1] is the skewness parameter and a = cos^. In this paper we 
consider positively skewed processes {9 = 1) say H^{t), i > 0, whose characteristic 
function writes 

^e-*^-*^)"' (1.35) 

where we used the fact that |^| e"*^'''^"^^) = -i^. The process H'^it), i > 0, has the 
important property of having non-negative, stationary and independent increments, 
and thus it is suitable to play the role of a random time. The law h^{x,t), x > 0, 
t > of H'^{t), t > 0, with Fourier transform h^{£^,t) and Laplace transform 

h,{fi,t) = e-*^", (1.36) 
solves the fractional diffusion equation, for ly E (0, 1], 

' im + -&) hA^^i) = 0' x>o,t>o, 

K{x,0) = S{x), (1.37) 
_K{Q,t) = 0, 

where the fractional derivatives are intended in the Riemann-Liouville sense. We 
notice that the process given by the composition of r S N independent stable 
subordinators of the same order i/ G (0, 1), say iH" (2^" (• • > H'^{t) ■ ■ •)), t > 
has law which reads 

Prjiff'- (2g'^(--vg'^(0---)) £ dx} ^ 

dx 

nOO nOO poo poo 

= / dsiih^{x,Si) / C;S2 2/li/(si, S2) / dS3 3h^{s2,S3) ■ ■ ■ dSr rhu{Sr,t). 

Jo Jo Jo Jo 

(1.38) 

In view of (1.35) and (1.36) we can easily write the Laplace and Fourier trans- 
forms of (1.38). For example the Laplace transform reads 

POO 

Ee-t^^H^2H^-^H-'(t)-)) ^ / dxe-'^^'PiiiH" {■■■rH''{t)---)) edx} 

Jo 

/•OO /"OO /"OO 

dsie"''!^ / ds22hu{si,S2) 3^I/(S2, S3) ' ' ' / dSrrK{Sr,t) 

Jo Jo Jo 

ds2e-'^''" / ds3 3K{s2,S3)--- dsrrhAsr,t) = e-*"'' (1.39) 



and therefore we have the following Fourier transform 



STABLE PROCESSES AT RANDOM TIME 



7 



1.2.2. Multidimensional stable laws and fractional Laplacian. Let us consider the 
isotropic n-dimensional process S'^^{t), t > , /3 G (0, 1], with density 

vp{x,t) ^ [ d^e-'«•^e-*"«"'^ a;eM",i>0, (1.41) 



(27r)" 

and therefore characteristic function 

dpitt) - Ee'«-«"'W = e-*ll€ll'', (1.42) 

where the symbol ||-|| stands for the usual Euclidean norm. The law (1.41) is the 
solution to the fractional Cauchy problem, for (3 e (0, 1] 

U§-^+{~Af)vp{x,t) ^ 0, xeW\t>0, ^^^^^ 
\vi3{x,0) = S{x). 

The fractional negative Laplacian appearing in (1.43) has been considered by many 
authors (see for example Balakrishnan [2]; Bochner [4]). The Bochner representa- 
tion of the fractional Laplacian reads 



- (-A)*^ = ^^^^^ / dAA^^-^A- A)"' A. (1.44) 
7r Jo 

Equivalently, an alternative useful definition can be given in the space of the Fourier 
transforms, as 

-i-Afu{x) ^ «||^||2^^i(0 (1.45) 

where 



Dom(-A)'' = |y e m ■■ j^^ \um' (l + IICII'^) di < (X)| 



(1.46) 



In the one-dimensional case and for < 2/3 < 1 we have that (see for example 
D'Ovidio ct al. [6] for details on this point), 

g2 Y ^2/3 

where g"^^2i3 is the Riesz operator usually defined as 

= ;t— ^w^^x: / Tzr^dz (1.48) 



d\x\^P 2cos/37rr(l-2/3)dx7_^ |a;-z|2/3 

and for which the Fourier transform becomes 



g2f! 



\C\''uiO- (1.49) 



2. Generalized fractional equations 

2.1. Linear combination of stable processes. In this section we start by con- 
sidering processes of the form 

m I 

:tf''i—'"(t) = ^Aj^fff (t), t>0,iy, e {0,1), j = !,■■■ ,m, (2.1) 

where Hj^{t), t > 0, are independent stable subordinators of order i^j g (0,1] 
introduced in section 1.2.1. Furthermore we will deal with the inverse process of 
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^vi,--- ,v,„ ^ say X''!' " ''^'"(t), t > 0, which can be defined as the hitting time of 
^uu -.u^ as 



(2.3) 



(2.4) 



The definition (2.2) of the process L'^^'"' '^"^ permits us to write 

PT{L''^ -^'''-{t) <x} = PT{^''^- -'''-ix)>t}. 
We present the following two results. 

Theorem 2.1. We have that 

i) The solution to the problem for Vj G (0, 1), j = 1, • • • , m, 

h^,^... ^^^{Q,t) = 0. 

is given by the law of the process " ^""(t)^ t > 0, defined in (2.1) 

ii) The solution to the problem for Uj e (0, 1), j = 1, • • • , m, 

is given by the law of the process L'^^'"' ''^"^(t), t> 0, defined in (2.2). 

The fractional derivatives appearing in (2.4) and (2.5) mwsi 6e intended in the 
Riemann-Liouville sense. 

Proof of i). Since for the Riemann-Liouville fractional derivative we have that 



(2.5) 



dx"" 



u{x) 



iO = HO'uix) 



(2.6) 



we can write the Fourier transform of the problem (2.4) as 



dt 



and therefore we have that 



dx 



(0 = E^^ HC)"^"^.!,...,.™ (C,*), 



V,-..™(C,0) = 1. 



(2.7) 



The Fourier transform of the law ky^^... ^^^^{x^i) of the process (2.1) is written as 

1 



for which 



7=1 A, (-ic)-^- 



(2.9) 



This is tantamount to saying that the Fourier transform of h^^_... ^^^^(x^t^ is the 
solution to the problem (2.7) and thus h^j^ ... ,j^{x,t) is the solution to (2.4). □ 
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Proof of ii). In this proof we will make use of the Laplace transform of the Riemann- 
Liouville fractional derivative which, in view of (2.6), can be written as 



u{t) 



(m) = /^""(m) ■ 



Taking the Laplace transform of (2.5) with respect to t we get 
™ ~ d ~ 

and performing the a;-Laplace transform of (2.11) we arrive at 

m 

The boundary condition appearing in (2.12) can be derived from (2.5) as 

poo , i^^j 

4.x.-..™ (0,M) = / d^e-^*^A, ^ = 5]A,m''^-i 



and thus from (2.12) we have that 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



Now we show that the Fourier-Laplace transform of the law of the process L'^^'"' (t), 
t > 0, coincides with (2.14). By taking into account the property (2.3) of the law 
of L^^'"' we obtain 

/•oo poo 



1 

M Jo 
1 

M Jo 

(1.36) 1 

M Jo 



dxe-'"" 
d 



dx f 



dx e 



dx 

' d_ 
dx 

d_ 
dx 



dx Jo 



dz fi^^^... ^i,^{z,x) 



(2.15) 



which coincides with (2.14). The proof of Theorem 2.1 is thus concluded. 



□ 



2.2. Generalized fractional telegraph-type equations. In this section we study 
equations of the form 

m C fiu 

(2.16) 



for Vj e (0, 1], j = 1, 



I, /3 e (0,1], c > 0, A > 0. The symbol 



stands 



for the Dzerbayshan-Caputo fractional derivative. Equation (2.16) generalizes the 
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telegraph equation in that an arbitrary number m of time-fractional derivatives 
appears and the n-dimensional fractional Laplacian governs the space fluctuations. 
Concerning the equation (2.16) we present the following result. 



Theorem 2.2. The solution to the problem for Vj G (0, 1], j = 1, • • • , m, /3 £ (0, 1], 

\wti,-- ,Urr,{x,Q) = 5{X). 

(2.17) 

is given by the law of the process 

W-i'-.--(f) = 5f (c2x''i'-''^"(i)) , t>Q, (2.18) 

where Sff is the isotropic vector process dealt with in section 1.2.2 and L'^'^'"' ''^'"{t), 
t > is the process defined in (2.2). 

Proof. Since for the Dzerbayshan-Caputo fractional derivative we have that, 



dt^ 



u{t) 



(m) - M''"(M)-M""'^i(0), i^e(0,l), 



(2.19) 



we can write the Laplace transform of (2.17) as 



Y.X,^i''^w'l^,..^,^{x,^,)-Y,\,^,'^^-H{x) = -c^ (-A)'^ ^ (a;, ^). (2.20) 

i=i j=i 
The Fourier-Laplace transform of (2.17) is therefore written as 



2 ||d||2/5' 



(2.21) 



Considering (2.3) we can derive the Fourier-Laplace transform of the process (2.18). 
We have that 



1/3 



P pGO poo 

(^,m) = / da;e'^'^ / dte^^^ j dsv/s {x,c'^s) C^,^... {s,t) 

JR" Jo Jo 







1 

M Jo 
M Jo 

(1^6) _ 1 /■ 

M Jo 



2 II II 2/3 

dse-'^ 



die"''* 
d 



d 



ds Jo 



^j,... ,„,„(z, s) dz 



2 II ^\\2B 

dse-'" 



ds 
d_ 
ds 



(2.21). 



(2.22) 



Since the Fourier-Laplace transform of the problem (2.17) coincides with that of 
the law of the process S-^^ (^c^L"^'"' '"""{t)), t > 0, the proof is complete. □ 
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2.3. Telegraph-type equations with two time-fractional derivatives. When 
in the equation (2.16) only two time derivatives appear we can rewrite the problem, 
for a,v € (0, 1] as 

r + 2A^) wl^ix, t) = -c2 (-A)'^ w^^ix, t), x€R-,t>0, 
\wi,{x,0) = d{x). 

(2.23) 

For a — 2v, v G (O, ^] the reader can recongnize in (2.23) the standard form of 
the classical fractional telegraph equation, investigated from a probabilistic point 
of view in Orsingher and Beghin [12] (for n = 1 and /? = 1) and in D'Ovidio et al. 
[6] (for n g N and /3 £ (0, 1)). In view of Theorem 2.2 is it not difficult to prove 
the following result. 

Corollary 2.3. The solution of the fractional Cauchy problem (2.23) is given by 
the law of the process 

W^'^it) = Sf {c^ L""-" {t)) , t > 0. (2.24) 

where 

L^'^^'it) = inf {s > : :?/"'''(s) = + (2A)^ H^{s) > , (2.25) 
for Hf and H2 independent stable subordinators. 

Proof. The proof of this result can be carried out by repeating the arguments of 
Theorem 2.2 and will not be reported here. It is sufficient to assume that Ai = 1, 
A2 = 2A, A > and A^ = for j > 2. □ 

2.4. The case a = kv. Let us consider a ^ kv, v ^ k ^ N, in (2.23). The 
problem becomes 

r(^+2A^)<,(a:,t) = -^{-/\)^wl^^{x,t), =r£M",i>0, 
UL..(^>0) = 5{x). 

(2.26) 

In view of Corollary 2.3 the solution to (2.26) is given by the law of the process 
S"^^ [c^L^^it)), t > Q. The Fourier-Laplace transform of w^^^{x,t) can be now 
written as 



k ,,1 k 



(2.27) 

where are the roots of ^"^ + 2Am'' + ||^ H^*^ = 0. 
For A: = 3 we get 

where A, B and C are the solutions to + 2Xfi'' + c^ = 0. Formula (2.28) 

can be rewritten as 

_ (m^"-! + 2A^'^-i) [/ 1 1 \ 1 
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1 



1 



1 



1 



^■^ - A fi^ - B J {A- B){B -C) 
1 



fi'' -A ^1'' -C J {A-C){B -C) 

1 



1 



1 



fi-" ~C {A~C){B ^C) 
By considering now the relationship 



^1'' - A{B - A){C - A) fi^ - B (A- B){C - B) 

(2.29) 



i/-(l-2iy) 



(2.30) 



we can invert (2.29) with respect to fj,. Thus we can exphcitely write the charac- 
teristic function of the process S^'^ [c^L^'^-''^{t)) , i > 0, as 



w: 



t-^^E.s^^u {At'') + 2Ag,,i [At'') (gt") + 2A£;,a (^t") 

(B - A) (C - A) (^ - B) (C ~ B) 

{A-C){B-C) 



(2.31) 



3. Multidimensional Gauss-Laplace distributions and infinite 

compositions 

In Orsingher and Beghin [13] the authors have shown that the process 

3„(i) = Bi (1^2 (IB3 • • • {\Bn+i[t)\) • • -1)1) , t > 0, (3.1) 

converges in distribution for n ^ 00 to a Gauss-Laplace (or bilateral exponential) 
random variable independent from < > 0. In this section we show that the process 
Bn (£^^' " (t)), i > 0, converges in distribution, for r — >■ (X), to a multidimensional 
version of the Gauss-Laplace r.v. and solves the equation, for Vj G (0, 1), r G N, 

Y.^3^^ti:-,ujx^i) = c2Am^;:....,^(a;,i), x e W,t> 0. (3.2) 
,=1 5* ^ 

The process ii^'' " '''"(^)i t > 0, is defined as 

= inf{s>0:jo^i'-''^'"(s)>t}, t>0 (3.3) 

where 

m 1 

^"''■■■'""W = 5] Aj^iiJ''^(2i/''^(3i?''M---ri?'^^W- ••))), i>0. (3.4) 

We start by presenting the following results. 
Corollary 3.1. We have that 
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i) The solution to the problem for Vj G (0, 1), j = 1, • • • ,m, r € N, 

K,,-.uJ^,0) ^ 5{x), " (3.5) 
is given by the law of the process 

in 1 

^J^'-'^-W - ^AfiiJ^^^CsH-'M 3ff"M--- r^^^^W- ■■))), *>0. (3.6) 

ii) The solution to the problem for Vj G (0, 1), j = 1, • • • , m, r g N, 



(3.7) 



is given by the law of the process 

^ > (s) > i}, i>0. (3.8) 

Proof of i). The proof is carried out in the same spirit of Theorem 2.1, thus by 
considering the Fourier transform of (3.5) we get 

The proof is completed by observing that the solution to (3.9) is given by the Fourier 
transform of the law of the process Jo^^' " ''^'"{t), t > 0, which can be obtained by 
means of the calculation 

1 

(3.10) 

that is the solution to (3.9). □ 
Proof of a). By considering the double Laplace transform of (3.7) we have that 

m 

E ^^^"^^"^I^i,- (^' A^) = (0' A^) - (^' A*) ' (3-11) 

where the boundary condition is given by 

m m 

i=i ^ J' i=i 
and thus 

(7,/i) = ' ^ (3.13) 

The definition (3.8) permits us to state that the processes S,'^^'"' -""^{t), t > 0, and 
^1^1, ••• ,1^™ (i)^ i > 0, are related by the fact that 

Pr ■ (i) < 4 = Pr{555:i'- '''"(x) >t}, (3.14) 
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and thus we can perform manipulations similar to those of Theorem 2.1. We have 
that the double Laplace transform of the law [^^ ... j^^{x,t) is then given by 

d 







dzKi,-,,.,„(.z,x) 



1 

M Jo 
1 

M Jo 



(lAO) 1 

M Jo 

and coincides with (3.13). 



dx f 



dx e'^"" 



dx t 



dx 

d_ 
dx 



Ee-M T.T=i ' iH"' ( ( 3H"3 (••• (t)-))) 



i=l ^3/^ ^ 

,=1 AjM ^ + 7 



(3.15) 

□ 



Theorem 3.2. T/ie solution to the problem for Vj G (0, 1], /3 G (0, 1] , = 1, • • • , 

/Er=iA,^roe;:...,.„(a3,i) = -c^ (-A)^ ,„(a;,t), a;eM",i>0, 

\^yC-,Vr„ (a;,0) = (5(a;) , 

(3.16) 

is given by the law 0/ the process 

21J„ {t) = Sff (c2£5;i'- ■''-(0) , t>0. (3.17) 

where the process fij^^' " (i), i > 0, is defined in (3.8). For (3 — 1, the process 
(3.17) becomes the subordinated Brownian motion B„ {c?Z':(}'"' '"^'"(i)), i > 0. 

Proof. The Fourier-Laplace transform of (3.16) can be easily derived as in Theorem 
2.2 and reads 



"'i.i,...,i.„.(^'Ai) 



By considering the law of the process S"^^ (c^H^^' " '"'"(t)) we have that 
dxe'^'' dte""' dsvp{x,c^s) il^ ... „^{s,t) 



(3.18) 







d_ 
ds 



_ 1 
M Jo 
1 

M Jo 

(1J39) 1 

M Jo 



213 f d 



' ^ Ee"^^7=i-^7' iff "^(2-ff''^(3fl"'^(-.H"^ (*)•••))) 
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ds e 



\2I3 



which coincides with (3.18). 



(3.19) 

□ 



We now consider the hmiting case for r — > cxd where the iteration of the process 
Sjf (c^£5^i (t)) , i > 0, is infinitely extended. In the next theorem we have 
that the limiting law of 



2P / 2ril/i,--- ,l^„ 



lim Si^ ic'£ 

r— J-oo 



t > 0, 



(3.20) 



is, for P — 1, a generalization to M" of the Gauss-Laplace probability density. This 
result represents an extension to the n-dimensional case of the infinitely iterated 
Brownian motion. 



Theorem 3.3. The distribution of the limiting process 

lim B„ (c^£-— >(t)) 

does not depend on t and reads 

Pr {X„i^n e dx} 



(-^1 5 ■ ■ ■ ; -^m) 

1 



(27r) 



dx 



\X\\ Kr. 



E"rA, 



The density (3.22) solves the equation 



92 



vi=i / i=i 



which is obtained from (3.16) by letting r — > 00. 
Proof. By assuming 



1 



the density 



Xo, 



B = 



3=1 



has first-order derivative which reads 



dx 



-tVmix) 



A 



n~2 



+ E- 



Kr. 



Bixl 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



„2\ 2 
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Kz 



„2 \ 2 



In the last step we applied the relationship 

of Lebedev [11], page 110. The second-order derivative now becomes 



2\2 



ABx]\Y^x] 



K 



AB'^x] 



K 



(b{x\ 

',{b{xI- 



2 



2 



2 



By considering the relationship 

K,+^{z) ^ K,^i{z) + 2-K,{z) 
z 

of [11], page 110, the derivative (3.28) takes the form 



2\ 2 



rom (a;i. 



K^,-,{B{x\ 



The Laplacian of rDm(a^i, • ■ ■ , a^m) therefore becomes 



J = l J 



/ 



{x)^AB^\Y^: 



Kr. 



and thus taking A and i? explicitely we obtain the desired result 



1 



I ('^l 1 ' ' ' : ^n) — ^ ^ -^j (-^1 ; ' ' ' ; ) 

i=i 
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□ 

Remark 3.4. For r — >■ oo the Fourier-Laplace transform (3.18) becomes 

1 y"-^ X, 

m) = ' ^, (3.33) 

^Er=i^.+c^ii^ii 

and thus the Fourier transform takes the form 

l^m(^) = ' ^. (3.34) 



The inversion of the Fourier transform (3.34) can be carried out by means of the 
hyperspherical coordinates. Thus we have that 

.0 = ^_ / ^^-1 



1 



(27r)" 

27r 



/ d6ii / (162- ■■ I dOn- 



^—ip[^l sin ^1 sin ■ " sin 9n-2 sin 0+a:2 sin 0i sin 02 ■ " sin 9n-2 cos ^ 
^— ip[2;3 sin 6i sin ^?2' " sin 6^ — ^ cos ^^ti-2+ +2^71-1 sin f?i cos cos ^i] 

sin"-2 6I1 sin"-3 02 ■■■ sin 6i„_2 ^^"^ ^ 



p^-^dp dOi d92--- d6i„_2 sin"-2 5)^sin"-3 5)2...sin6'„_2 



(27r)"-i 

^—ip[x^ sin ^1 sin ^2' " sin ^n,_3 cos 0^1-2+ -\-^n — i sin cos cos ^i] 

E™iA, 



Jo pA/.Tf + a;^ sin6'i sin6'2 • • • sin6'„_2 (3.35) 
Lj=i Ai + V ^ / 

We now evaluate the integrals with respect to 9j by means of formula 6.688 page 
727 of Gradshteyn and Ryzhik [9], which reads 

/ sin'^^"^ a; cos (/Scosx) Jjy (asinx) da; = W — - — r>^i/+i f\/a^~+^V 

Jo V 2 (a2 + ^2)2 + 4 2 V J 

(3.36) 

valid for > — 1. We start with the integral with respect to 6'„_2 

J (i6»„_2 e-*P^3 sin Si . ■ . sin e„_3 COS e„_2 sin e„_2 Jo ^p^x^ + ^2 sin 6*1 ■■ ■ sin 6'„_2 j 

= 2 y d6'„_2 cos (pzs sin 6*1 ■ ■ ■ sin^n-a cos6'„_2) sin6'„_2-'o (^P^J ^\ + sin 6*1 ■ ■ ■ sin6'„_2^ 



= v'Stt (^psin^i • ■ • sinS,i_3y + + XgJ Ji |^Py + + ^3 sin^i ■ ■ ■ sin9„_3 

and thus the integral with respect to 0„_3 becomes 
I dOn-se-'P'''''''''''-'''''^"-'''"''''"-' 



psin^i • • • sin0„_3Wa;2 + a;2 + ajg j Ji P\/ + a;2 + a;^ sin^i • • • sin 



7)1-3 
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2v27r I psin0„_i • • • sin6'„_4W2:^ + 2:2 + I / c?0„_3 sin^ 0„_3 



cos {px4 sin 6*1 • • • sin6'„_4 cos On-3) Ji yP\J x\+ x\+ a;| sin 6*1 •• • sin0„_3 

/ \ 2 _ 1 

= f V2tt\ {p^ sit? 6*1 •• • sin^ 6'„_4 {xl +xl+xl+x'^^) " 

Ji (^\J x\ + x\ + x\ + x\ sin 9i sin 02 ■ • • sin 6'„_4^ . (3.37) 
After n ~ 2 integrations we arrive at the integral with respect to p which reads 



tvi {x) = (27r) 



E 



x2 



which, for (3 = 1 and after the change of variable p^c^ = y^, becomes 



(3.38) 



(27r)-ti 
c 



E;.r-y 



for n<5 



(27r)t 



(27r) 



e;LiA, 



e;1iA, 



EjlrA, 



|a;|| K. 



E,"lrA, 



(3.39) 



where we used formula 6.566 page 679 of Gradshteyn and Ryzhik [9], which reads 

/ 

Jo 



1 3 

dxx^+^JJax)^ = WK^(ab), a > 0, 5R(6) > 0, -1< •Riv) < -. 

+ 0-^ 2 



(3.40) 

Remark 3.5. We can check that (3.22) for all n e N is a true probability density. 
Wm{x)dx = 



area 



(Sf) / \J^T=i 



dp 



(27r)t 1 / yEjliA 



r f (27r). 



P^KlL^l p 



dp = 1 (3.41) 



in force of formula 6.561(16) of Gradshteyn and Ryzhik [9] page 676 



x^^K^ (ax) = 2^'-^a-''-^T 



(3.42) 
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valid for 3? + 1 ± z^) > and 3?(a) > 0. The non-negativity of (3.22) is shown by 
the following integral representation 

/>oo 

K^{z) = / e-"™''''*coshi/tdi (3.43) 

valid for |arg(z)| < f (see [9] page 917 formula 8.432). 
By considering that 

K_.^{z) = Kr{z) = y^e-^ (3.44) 
from (3.22) we derive the following probability density for a; e M'^, 

(2c)^^E;.i-f {2c)H\\xf ^'-"'^ 

In the two dimensional case the distribution (3.22) has a simple structure which 
reads 

n,mixi,X2) - ^ ' . (3.46) 

In view of (3.44) it is also easy to show that the distribution (3.22) coincides for 
n = 1 with the classical Gauss-Laplace distribution. We have that for n = 1 (3.22) 
becomes 

3 




l"l (3.47) 



2c 

Furthermore, for Ai = 1, A2 = 2A, A > and Xj = for j = 3, • • • , m, we note that 
(3.47) coincides with formula (3.18) of Orsingher and Beghin [13]. 

Remark 3.6. By considering the iterated random walk 

r„(fc) = 5i (52 (•••(^„(fc)) •••)), fceN, (3.48) 

with Sj, j — I, ■ ■ ■ ,n, independent random walks, Turban [17] has shown that for 
n — >■ 00, Yn{k) converges to a stationary r.v. (independent from k) which possesses 
Gauss-Laplace distribution, in accord with result (3.22) of the present work and 
with (3.12) of Orsingher and Beghin [13]. 
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